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Abstract 

In this talk, I present my personal view on the status 
of lattice QCD calculations. I emphasize the role played 
by the chiral perturbation theory (^PT) in analyzing the 
lattice data of various physical quantities, including chi- 
ral condensate, topological susceptibility, and pion mass 
and decay constant. I then discuss on the status of de- 
termination of fundamental parameters and other quan- 
tities of phcnomenological interest. 

1. Introduction 

It was already more than 30 years ago that Quan- 
tum Chromodynamics (QCD) was proposed as the fun- 
damental theory of strong interaction. Since then, there 
have been an enormous number of experimental data 
that support QCD including its quantum effects. Those 
experiments are essentially measuring pcrturbative as- 
pects of QCD at high energy, while the non-pcrturbative 
dynamics at low energy still remains as a difficult prob- 
lem. Solving QCD is difficult mainly because its vac- 
uum is so complicated. There is no single dominant 
gauge field configuration; it is not completely random 
cither. Since hadrons arc floating on this vacuum, an 
understanding of the QCD vacuum is a prerequisite to 
calculate hadron properties starting from first principles. 

Lattice QCD is one of the methods to regularize ultra- 
violet divergences in QCD. Unlike the commonly used di- 
mensional regularization, whose definition involves per- 
turbation theory, lattice QCD is valid in both ultraviolet 
and infrared regimes. As the regulated theory is math- 
ematically well-defined, direct numerical calculation of 
the path integral that defines QCD is possible. Thus, lat- 
tice QCD provides a powerful method for first-principles 
calculation of QCD including its non-pcrturbative dy- 
namics at low energy. It requires huge computational 
resources, especially to incorporate the fermion loop ef- 
fects in the path integral. Since 1980s, lattice QCD has 
been using the high-end supercomputers available at the 
time. It is remarkable that researchers of lattice QCD 
have even developed machines that lead the entire su- 
percomputing. 

The state-of-the-art of lattice QCD simulations could 
be summarized as follows: "Realistic simulation of QCD 
to study the static properties of hadrons is now feasible" . 
It means that the inclusion of up, down and strange 
sea quarks has become a standard at small enough lat- 
tice spacing (a < 0.1 fin) on large enough volume (L > 
2.5 fm) to hold a single hadron. The low-energy hadron 
spectrum, for instance, has been well reproduced by sev- 
eral lattice groups. 

Although the recent progress of lattice QCD is so im- 
pressive, it still has many limitations. One of those is the 
multi-scale problems. The scales that may enter in the 
QCD phenomena are widely ranged: up and down quark 
masses (~ 5 McV), strange quark mass (~ 100 MeV), 
the QCD scale (~ 300 MeV), charm quark mass (~ 



1.5 GeV), and bottom quark mass (~ 4.5 GeV). Treat- 
ing light quarks requires more computational costs that 
grows as 1/m 2-3 with m q the light quark mass; reduc- 
ing the lattice spacing to treat heavy quarks needs more 
resources that typically scales as 1/a 7 with a the lat- 
tice spacing. Therefore, the best strategy for practical 
applications is to use "effective theories" such as the chi- 
ral perturbation theory (%PT) for light quarks and the 
heavy quark effective theory (HQET) for heavy quarks. 
The light and heavy quark masses for which these effec- 
tive theories are valid have to be carefully investigated 
using lattice QCD calculations. 

This talk is not a comprehensive review of the field. 
Rather, I would like to present my personal view of the 
status of lattice QCD. I start the discussion from a study 
of the fundamental property of the QCD vacuum in Sec- 
tion 2. Then I summarize the calculations of some inter- 
esting phcnomenological quantities of light hadrons with 
special emphasis on the convergence property of the chi- 
ral expansion (Section 3). Determination of the funda- 
mental parameters of QCD, such as the strong coupling 
constant and quark masses is of particular importance, 
that I discuss in Section 4. My discussion on heavy flavor 
physics is brief (Section 5). My summary and perspec- 
tive are in the last section. 



2. Chiral symmetry breaking 

2.1. Chiral symmetry and lattice QCD 

One of the fundamental properties on the QCD vac- 
uum is the spontaneous breaking of chiral symmetry. 
Even before QCD, many important properties of low- 
energy hadrons, such as the GMOR relation and other 
soft pion theorems, were discovered based on the PCAC 
relation and current algebra. From the modern perspec- 
tive, they arc derived from the chiral effective theory, 
which is constructed assuming the spontaneous chiral 
symmetry breaking. For the thorough understanding of 
strong interaction, it is therefore a crucial step to estab- 
lish a link between QCD and chiral effective theory. 

Chiral symmetry of course plays a key role in the un- 
derstanding of chiral symmetry breaking. In the fla- 
vor non-singlet sector of chiral symmetry, pions arise as 
the Nambu-Goldstone boson associated with the sponta- 
neous symmetry breaking, while in the flavor-singlet sec- 
tor the chiral symmetry is violated by the axial anomaly 
and is related to the topology of non-Abelian gauge the- 
ory. There are near-zero modes of quarks associated 
with the topological excitations; their accumulation in 
the vacuum leads to the symmetry breaking in the fla- 
vor non-singlet sector as indicated by the Banks-Cashcr 
relation pQ. Therefore, the initial setup to study the chi- 
ral symmetry breaking should preserve the flavor singlet 
and non-singlet chiral symmetries. 

There is a problem in realizing the chiral symmetry 
on the lattice. The conventional Wilson-type fcrmions 
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violate the chiral symmetry at the action level, and the 
discrimination between the physical effect of symmetry 
breaking and the lattice artifact is not clear. On the 
other hand, the staggered fermions have a chiral sym- 
metry but break the flavor symmetry. With these lat- 
tice fermions, the continuum limit has to be taken before 
analyzing the data using the continuum chiral effective 
theory. 

The domain- wall [2J EJ S] and overlap [5] [B] fermions 
solve this problem in a theoretically clean manner. Since 
they satisfy a modified version of chiral symmetry at fi- 
nite lattice spacing [TJ, the continuum-like axial- Ward- 
Takahashi identities are hold on the lattice. (The chi- 
ral symmetry is exact for the overlap fermion; for the 
domain-wall fermion the chiral symmetry is restored in 
the limit of infinite size of the fifth dimension, but small 
violation remains in practical simulations.) This means 
that the soft-pion theorems are all satisfied on the lat- 
tice just as in the continuum theory, as far as the chiral 
symmetry is spontaneously broken. The use of X-PT is 
therefore justified at any finite lattice spacing a. (For 
the Wilson and staggered fermions, one has to include 
the terms that represent the violation of chiral symme- 
try) Although the numerical cost is substantially high 
compared to the Wilson or staggered fermions, these 
fermion formulation should therefore be used when the 
chiral symmetry is crucial. 

The spontaneous chiral symmetry breaking is probed 
by the chiral condensate (qq). Its lattice calculation is 
challenging because the scalar density operator qq has 
a power divergence of the form m q /a 2 as the cutoff 1/a 
goes to infinity. The massless limit has to be taken to ob- 
tain physical result. (When the chiral symmetry is vio- 
lated from the outset as in the Wilson-type fermions, the 
divergence is even stronger ~ 1/a 3 -) On the other hand, 
the condensate vanishes in the massless limit, when the 
space-time volume is kept finite. Therefore, the proper 
order of the limits is to take the infinite volume limit first 
and then the massless limit, which is called the thermo- 
dynamical limit. 
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Fig. 1 Spectral density of the Dirac operator calculated by the 
JLQCD collaboration (histogram). The curve represents a 
next-to- leading order xPT calculation 1111 . The pion mass is 
around 300 MeV and the system in the p-regime. 
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Fig. 2 Chiral condensate as a function of sea quark mass. Re- 
sults from the JLQCD collaboration. 



2.2. Spectral density and chiral condensate 

The problem of the ultraviolet divergence can be 
avoided by focusing on low-lying eigenmode spectrum of 
the Dirac operator. As indicated by the Banks-Casher 
relation [T], the chiral symmetry breaking is induced 
by an accumulation of low-lying eigenstates of quark- 
antiquark pair 
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where p(X) denotes the eigenvalue density of the Dirac 
operator, p(\) = (1/V^) ^ fc (<5(A — A^)}. The expectation 
value (• • • ) represents an ensemble average and k labels 
the eigenvalues of the Dirac operator on a given gauge 
field background. On the right hand side of QJ, £ is the 
chiral condensate, £ = —(qq), evaluated in the mass- 
less quark limit. In the free theory, we expect a scaling 
p(X) ~ A 3 for a dimensional reason and thus p(0) = 0. 
The relation fl} implies that the spontaneous chiral sym- 
metry breaking characterized by non-zero S is related to 
the number of near-zero modes in a given volume. 

Based on xPT, more detailed forms of p(X) at finite 
A, V and m have been obtained. This provides a theo- 
retical basis to control the scaling under the change of 
these parameters. A recent simulation of the JLQCD 
collaboration gave the spectral density with 2+1 flavors 
of dynamical overlap fermions [8]. (An earlier analysis 



is in [HI [IH])- One of the results is shown in Figure [I] 
which is obtained at lattice spacing ~ 0.11 fm and lat- 
tice volume ~ (1.8 fm) 3 x(5.4 fm). The plot compares 
the result with a x?T calculation at the next-to-leading 
order [llj . The shape of the spectrum is nicely repro- 
duced, and the height determines the chiral condensate 
at the given quark mass. 

A similar analysis is attempted using the Wilson 
fermion in [12]. Due to the explicit violation of chi- 
ral symmetry, the spectrum of the near-zero modes is 
distorted. Relatively larger eigenvalues are, therefore, 
mainly used in the comparison with x?T. 

An extrapolation of the chiral condensate Y*(m u d, m s ) 
to the chiral limit of up and down quark masses is shown 
in Figurc[2] The lattice data of the JLQCD collaboration 
show a curvature, wh ich is essentially a pion- loop effect 
as predicted by x?T [T3] 
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The data point close to the chiral limit, which is in the 
so-called e-regime, is helpful to identify this curvature. 
A preliminary result in the chiral limit of up and down 

quarks is E m (0,m s ; 2 GeV) = [243(4)(+^) MeV] 3 . 
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Fig. 3 Topological susceptibility \t as a function of sea quark 
mass. Lattice data are from I15I I16I . 



2.3. Topological susceptibility 

The topology of gauge field configuration plays a key 
role in the accumulation of the near-zero modes. The 
amount of topological excitations is measured by the 
topological susceptibility xt = (Q 2 )/V, where Q is the 
global topological charge in a given volume V. An expec- 
tation from xPT is that xt behaves as \t = mS/7V/ with 
m the quark mass of Nf degenerate flavors [H]. (There 
is also a formula for non-degenerate quark masses.) A 
lattice result by the JLQCD and TWQCD collaboration 
[El [16] obtained with the overlap fermion is shown in 
Figure [31 The results from two- flavor and 2+Tflavor 
QCD are well described by x?T. This precise calcula- 
tion was made possible by a new method to calculate 
Xt through a topological charge density correlation on 
gauge configurations at fixed global topological charge 

Through these studies, the spontaneous breaking of 
chiral symmetry is well established using the first- 
principles calculation of lattice QCD. The exact chiral 
symmetry provided by the overlap fermion played a cru- 
cial role there. The x?T is confirmed to be valid near the 
chiral limit for fundamental quantities such as the chi- 
ral condensate and topological susceptibility. The next 
question would be how much the region of xPT is ex- 
tended towards wider applications and larger values of 
pion masses/momenta. 

3. Light hadron phenomenology 

3.1. Convergence of chiral expansion 

At Lattice 2002, the annual conference on lattice field 
theory, there was a panel discussi on o n the issue of chi- 
ral extrapolation of lattice data [18]. The problem at 
that time was that the curvature due to the chiral loga- 
rithm of the form m 2 In m 2 was not visible in lattice data 
for any physical quantity. This was mainly because the 
quark mass in the dynamical fermion simulations at that 
time was too large that pion mass was above 500 McV, 
which is presumably out of the range of xPT. This lead 
to a large systematic error in the chiral extrapolation. 

Since then, by the development of algorithms and ma- 
chines, the pion mass in the lattice simulations has been 
reduced down to 200-300 MeV. For a summary of recent 
large-scale simulations, see the plenary talk by Scholz at 
Lattice 2009 [T5]. Now, it is therefore the time to inves- 
tigate the convergence property of the chiral expansion. 

The xPT provides a systematic expansion in terms of 




Fig. 4 Comparison of chiral expansion in terms of x, x and £. 
The plots represent m^/m q (left) and (right). Fits of the 
three lightest data points with the NLO ChPT formulae l[3j 
and are shown. Results are from |20| . 



small to 2 and p 2 , but the region of convergence of this 
chiral expansion is not known a priori. With the exact 
chiral symmetry, the test is conceptually cleanest, since 
no additional terms to describe the violation of chiral 
symmetry has to be introduced. (With other fermion 
formulations, this is not the case as already noticed. The 
unknown correction terms are often simply ignored.) 

For the pion mass and decay constant / w the ex- 
pansion is given as 



= 2B 

Tlq 

U = f[l-xlnx 
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where and denote the quantities after the correc- 
tions while m and / are those at the leading order. The 
expansions ([3J and ([¥]) may be written in terms of either 
x = 2m 2 /(4^/) 2 , x = 2m 2 /(4^/) 2 , or £ = 2m 2 /(47rA) 2 
(we use a notation of f n = 131 MeV). They all give 
an equivalent description at this order, while the conver- 
gence behavior may depend on the expansion parameter. 

Among other works, I use our own data (by the 
JLQCD and TWQCD collaborations) for a discussion 
here. Figure [4] shows the comparison of different expan- 
sion parameters in two-flavor QCD [20]. The fit curves 
are obtained by fitting three lightest data points with 
the three expansion parameters, which provide equally 
precise description of the data in the region of the fit. 
If we look at the heavier quark mass region, however, 
it is clear that only the ^-expansion gives a reasonable 
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Fig. 5 fn/fn from various lattice groups. A plot from E. Scholz 
|19| . References for individual data can be found there. In 
the plot, PDG'08 is from the experimental data of leptonic 
kaon decay, assuming \V U s\ = 0.2255(19). 



function and others miss the data points largely. This 
clearly demonstrates that at least for these quantities 
the convergence of the chiral expansion is much better 
with the ^-parameter than with the other conventional 
choices. This is understood as an effect of resummation 
of the chiral expansion by the use of the "renormalized" 
quantities m\ and / w . In fact, only with the ^-expansion 
we could fit the data including the kaon mass region with 
the next-to-next-to-leading order (NNLO) formulae [20] . 

With the NLO formulae (like those in (J3J or (gj), the 
convergence of the chiral expansion is marginal in the 
kaon mass regime. In fact, some groups concluded not 
to use the SU(3) xPT for the kaon sector but use the 
SU(2) formula with the strange quark treated as a heavy 
particle. This corresponds to an expansion in terms 
of m u d/m s and is a theoretically consistent treatment, 
though the predictive power of xPT is lost to some ex- 
tent. A summary of the results for fx/ fir from a review 
talk at Lattice 2009 by Scholz [19) is shown in Figure [5] 
Different groups take different strategies on the treat- 
ment of strange quark in the chiral fit, i.e. SU(2) or 
SU(3), NLO or NNLO. At the level of the error of order 
±0.02, the results are in agreement. 



3.2. Kaon semi-leptonic form factor 

Whether the strange quark can be treated within the 
SU(3) is an important issue, since the main sources 
of phcnomcnological information for kaon physics relies 
on xPT. For instance, the determination of a CKM ma- 
trix element \V US \ uses the semileptonic decays K — ► irtu. 
It is made precise because the relevant form factor 
f + (q 2 ) is normalized to 1 at q 2 = in the degenerate 
SU(3) limit. The devi ation from there is estimated us- 
ing the SU(3) xPT [21]. Even the lattice calculations use 
its formula as a guide to fit the data. The recent result 
by the RBC-UKQCD collaboration with 2+1 flavors of 
domain-wall fermions [22] is / + (0) = 0.964(5), which is 
consistent with 0.961(8) of the early phenomenological 
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Fig. 6 Fit of the vacuum polarization function on the lat- 
tice with continuum perturbative calculation. Data from the 
JLQCD collaboration with 2+1 flavors of dynamical overlap 
fermions 1331. 
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3.3. Neutral kaon mixing 

Results for Bk, the bag parameter to characterize 
the K° — K° mixing, has been made very precise and 
stable by recent lattice calculations. The recent results 
in 2+1-flavor QCD are: B K = 0.720(13)(37) with 2+1 
flavors of domain- wall fermions [23], 0.724(8) (28) with 
domain- wall vale nce quarks over 2+1 flavors of stag- 
gered sea quarks [241. E ven including two- flav or calcula- 
tions using overlap [25| and twisted-mass [26] fermions, 
all the dynamical lattice calculations give consistent re- 
sults. The average given by Lubicz at Lattice 2009 [27] 

is B K = 0.731(7)(35). 

With the 5% error, Bk is no longer the dominant un- 
certainty in drawing the constraint on the CKM uni- 
tarity triangle from the indirect CP violation parameter 
\ck\- Since \zk\ is related to the (/?, rf) parameters of the 
CKM unitary triangle as |e#| oc B K A 2 i](c 1 +c 2 A 2 (l-p)) 
(ci and C2 are known constants), the uncertainty in 
\V c b \ = AX 2 is now more important, as emphasized in 

For the lattice calculations related to the kaon physics, 
I also refer the readers to a recent c omp rehensive review 
by Boyle at the Kaon09 conference [29] . 

4. Fundamental parameters 

Next, I discuss the recent determinations of fundamen- 
tal parameters in QCD, i.e. the strong coupling constant 
and quark masses. 



4-1- Strong coupling constant 

The determination of the strong coupling constant is 
equivalent to an input of the lattice scale 1/a. This is 
because, for a given lattice with a bare coupling con- 
stant a), at , it gives the running of the coupling constant 
CKg (a -1 ). But, in order to obtain the value in more fa- 
miliar definitions, such as the one in the MS scheme, one 
has to convert the result using perturbation theory as 

a^ ls (fi) = Z(^a)a' s at (a _1 ). Therefore, precise determi- 
nation requires a good control of systematic errors in the 
perturbative matching factor Z([ia). This is the reason 
that the previous lattice calculation by the HPQCD col- 
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Fig. 7 Summary of the quark mass determinations from the lattice. A plot from Scholz I19I . References for individual data can 
be found there. Results with perturbative and non-perturbative renormalization factors are distinguished by pluses and starts, 
respectively. 



laboration |30l I31] employed an automated perturbative 
calculation methods to calculate two-loop contributions 
on the lattice. 

Recently, new methods have been proposed. They are 
based on a physical quantity which is sufficiently short- 
distance and is calculable in both continuum and lattice 
theories. A well-known example is the Adler function, 
which is a derivative of the vacuum polarization function 
D(Q 2 ) = -Q 2 dU v (Q 2 )/dQ 2 . Since this quantity docs 
not have an ultraviolet divergence, the perturbative cal- 
culation in the continuum theory using the dimensional 
regularization can be applied for the lattice data without 
modifications except for possible discretization effects. 
Therefore, one may fit the lattice data using the contin- 
uum formula known to three-loops (or even four-loops for 
the leading term) supplemented by an operator product 
expansion in l/Q 2 . Then, one can directly determine 

aiF(At) PJ. The fit of lattice data in 2+1-flavor QCD 
is shown in Figure [6] for different light quark masses. A 
preliminary result a s (M z ) = O.llSl^X^Xlg) has an 
error which is compa rable to the previous best lattice 
calculation, 0.1183(8) (30] El]. 

A closely related method is the use of the charmonium 
two-point function. By taking an appropriate moment 
in the coordinate space, the ultraviolet divergence is re- 
moved and the lattice data can be directly fitted with 
a corresponding conti nuu m calculation. The result is 
a s {M z ) = 0.1174(12) [34]. With this method, one can 
also determine the charm quark mass at the same time 

M- 

Overall, the lattice calculation now provides the most 
precise determination of the strong coupling constant. 
Along the lines of the recent calculations, further im- 
provement is expected. 



4-2. Quark masses 

Light quark masses can be extracted with pion and 
kaon masses as inputs. Therefore, the results depend 
on the fit functions (SU(2) or SU(3), NLO or NNLO) 
and the mass range used in the analysis of pion and 



kaon masses. Since the pole mass of quarks is not well 
defined due to confinement, one has to use a renormal- 
ization scheme given at short distances such as the MS 
scheme. Conversion from the lattice bare quark mass is 
done using perturbation theory or partially using some 
non-perturbative methods. The errors with purely per- 
turbative method could be significantly underestimated, 
as the coefficients of higher order terms are not known 
a priori. The res ults summarized in Figure [7] (from the 
review by Scholz [19] at Lattice 2009) should be viewed 
with these caveats in mind. Since many results appeared 
only recently, the understanding of these systematics will 
be achieved in the coming years. I expect that the error 
is reduced to 2-5% in the near future. 

The isospin breaking, or the up and down quark mass 
ratio, cannot be determined solely from the QCD cal- 
culation, as the electromagnetic effect gives a substan- 
tial contribution. For instance, the charged and neutral 
mass difference of pion is dominated by the QED ef- 
fect, which has to be subtracted to determine m u — m e i- 
There is an attempt to identify the QED effect us- 
ing the Weinberg- type sum rule |35j : m 2 ± — m 2 cx 
olem J °° dQ 2 Q 2 [U v {Q 2 )-n A {Q 2 )]. Here, the difference 
of the vacuum polarization functions in the vector and 
axial-vector channels is involved, which means that this 
mass difference is triggered by the spontaneous chiral 
symmetry breaking. With exact chiral symmetry, lattice 
calculation of Uv{Q 2 )—Ha(Q 2 ) is possible [36]. Another 
interesting method to calculate the QED effect is a sim- 
ulation of the whole QCD+QED system [37] [38]. Then 
one can also calculate the charged and neutral kaon mass 
difference, which is useful to study the amount of the vio- 
lation of the Dashen's theorem rn 2 K± —m 2 K0 = m 2 . ± —m 2 2 
[39] . More striking application would be the mass dif- 
ference between proton and neutron, as it is related to 
a big question of why the ma tters in Nature are stable. 
For a recent analysis, see [38] . 
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Fig. 8 Inclusive and exclusive determinations of |V c j,| (top 
panel) and |V u j,| (bottom panel). Plats are taken from Van 
de Water [23]. The lattice data are those in I4UI . 



5. Heavy flavors 

The biggest motivation to consider the heavy quark 
flavors is to put constraints on the CKM unitarity trian- 
gle. Since the experiments of B meson decays have been 
dramatically improved over the last several years, the 
required precision for the lattice calculation has become 
high. Namely, to be interesting, the lattice calculation 
for the decay constants, bag parameters, and form fac- 
tors has to be typically as precise as 5% or even better, 
which is a challenge for the lattice QCD community. 

For instance, determination of the CKM matrix ele- 
ments \V c b\ and \V u b\ can be done through inclusive or 
exclusive decay modes. The inclusive determination uses 
pcrturbative calculation relying on the quark-hadron du- 
ality assumption, while the exclusive determination re- 
quires the lattice calculation of the form factors. So far, 
the precision of the inclusive determination is slightly 
better for both \V c b\ and \V u b\, as far as the error is taken 
as a face value. S ee a summary plot in Figure [5J made 
by Van de Water [55]. More importantly, the inclusive 
and exclusive determinations are inconsistent with each 
other at around 2a level. In order to understand this 
discrepancy, more precise calculations and experiments 
are required. 

Since the Compton wavelength of heavy quarks is too 
short to treat on the lattice that is available with present 
computational power, one has to use some effective the- 
ory. The idea is to factor out a trivial heavy quark mass 
dependence e~" nQt from physical amplitudes; remaining 
terms are organized as an expansion in terms of I/tuq. 
The effective theories realize this in a language of effec- 
tive lagrangians, and it is also possible to construct a lat- 
tice version. The problem is that the effective theory has 
many parameters that have to be determined so that it 
reproduces the original QCD. In many cases, the match- 
ing of these parameters is done by perturbation theory, 
with which neglected higher order corrections become a 
source of systematic error. Non-perturbative parameter 
matching has also be attempted |41j , but is not common 
as it required dedicated study for each choice of lattice 
gauge and fermion actions. (For an int roductory article 
on heavy quarks on the lattice, see [42]). 

On the other hand, as the available computational re- 
sources increase, it has also become possible to directly 
simulate heavy quarks with the lattice spacing kept as 
small as possible. It is feasible only in the charm quark 
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Fig. 9 ^(s) meson decay constant. Experimental data and 
lattice calculations are compared. A plot from Van de Water 
[281 . CLEO and FNAL/MILC data have been updated in 
2009; the solid lines arc the latest results. 



mass region, and an extrapolation is necessary toward 
bottom quark mass. In particula r a v ariant of the stag- 
gered fermion (the HISQ action) [43] has been found to 
be useful to minimize discretization effects. Because the 
staggered fermion has a chiral symmetry (at a price of 
broken flavor symmetry), the order of discretization ef- 
fect is limited to 0(a 2 ), 0(a 4 ), the HISQ action is 
designed to remove the leading 0(a 2 ) effect so that the 
remaining error is 0(a i ) at the tree level. Other non- 
chiral fermions have the errors of 0(a), 0(a 2 ), the 
removal of the lowest order effect leaves the Oia 2 ) ef- 
fect, which is still substantial. 

In order to demonstrate how precision is crucial in 
heavy flavor phenomenology, I show a summary plot of 
the calculation of D( s ) meson decay constant /n (s) in 
Figure [5] Although there are three recent lattice calcu- 
lations, the result of the HPQCD collaboration [44] using 
the HISQ action has an order of magnitude smaller error 
quoted and thus dominates the lattice calculations. If it 
is compared with experiments, one may conclude that 
there is a significant discrepancy in fjj g , which triggered 
theorists to consider some new physics models that may 
explain this [5S1[3S]. Therefore, it is very important to 
have other calculations that reach this level of precision 
with different lattice fermion formulations. This needs 
two or three more years in my opinion. 



6. Summary and perspective 

The development of supercomputers was initiated in 
1980s, and an exponential growth of the computational 
power has been continued since then. The rate is about 
a factor of 10 in 5 years over the last 25-30 years, and 
still no trend of speed-down is observed. This, of course, 
has driven the rapid improvement of lattice calculations. 
Until 1990s, only the quenched calculations were possi- 
ble, and the results were subject of uncontrolled system- 
atic errors. Large-scale dynamical fermion simulations 
were started at around late 90s, and with lots of efforts 
to improve algorithms and techniques the simulations 
in a realistic setup have become feasible recently. This 
means that the lattice calculation can now produce real 
predictions from QCD for many interesting quantities. 

Theoretically, the biggest achievement in the last 10- 
20 years in lattice field theory is the formulation of lat- 
tice fermions with chiral symmetry. It removed a large 
class of limitations of lattice calculations. Indeed, spon- 
taneous chiral symmetry breaking and related phenom- 
ena can now be directly simulated in lattice QCD. 

These facts promise a wider range of applications of 



lattice QCD. The coverage of this talk has been lim- 
ited, but lattice QCD may be useful in many other areas 
of particle and nuclear physics. They include nucleon 
structure, spin physics, exotic hadrons, muon g — 2, fla- 
vor (B, D, and K) physics, heavy ion physics, and even 
dark matter search or neutrino interactions — nearly all 
the subjects that were covered by the Physics in Colli- 
sion conference. Lattice QCD is your friend! 
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